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NATIONAL ADVISOEY COMMITTEE EOE AEEONAUTICS 
EESEAECH MMOEANDUM 
COOLING OE GAS TURBINES 
IV - C ALCULATED TEMi’ERATUEE DISTRIBUTION IN 
THE TRAILING PART OF A TURBINE BLADE 
USING DIRECT LIQUID COOLING 
By ¥. B3rron Brown and Nilliam E. Monroe 


SUMMARY 

A "theoreLical analysis was made "to determine the temperature 
distribution in the trailing portion of a liquid-cooled turbine blade 
between the coolant passage and the trailing edge and to study the 
influence of various design and operating variables on the hot-spot 
temperature at the trailing edge. The trailing portion of a typical 
turbine blade was approximated by equivalent rectangular and wedge- 
shaped sections to facilitate analysis, and representative values 
for the boundary temperatures and surface heat- transfer coefficients 
were assumed or computed according to standard methods. 

Three-dimensional temperature distributions were obtained at 
effective gas temperatures of 2000°, 3000°, and 5000° F for a 4-inch 
prismatic section, cooled over the entire length and having a thermal 
conductivity of 15 Btu per hour per foot per °F, and at 2000° F for 
similar sections having thermal conductivities of 120 and 210 Btu 
per hour per foot per °F. In addition, comparison of one-dimensional 
temperature distributions through the outer portion of the blade was 
made for two shapes, equivalent rectangular and wedge plane sections, 
at an effective gas temperature of 2000° F. The influence of thermal 
conductivity on one -dimensional temperature distribution in both 
plane sections was evaluated. With the equivalent rectangular solu- 
tion, an investigation was made to determine the effect on trailing- 
edge hot-spot temperature of variations in the trailing-section 
width, thickness, thermal conductivity, and coolant-passage heat- 
transfer coefficient. 

The three-dimensional and one-dimensional solutions of the 
equivalent rectangular trailing portion give almost identical temper- 
ature distribution in the blade section away from the influence of 



2 


MCA EM No. E7Blld 


rim cooling. For tile blades of low thermal conductivity and wide 
trailing sections the hot-spot temperature closely approaches the 
effective gas temperature, and infinite increase in the surface heat- 
transfer coefficient at the coolant passage relative to that obtained 
with water will have small effect on the trailing-edge temperature. 
Further reduction of trailing-edge hot-spot teraperature can be 
obtained by the following three methods given in order of decreasing 
effect; by use of materials having high thermal conductivity, by 
altering the blade shape to provide a relatively short trailing 
portion, by locating additional cooling passages a short distance 
from the trailing edge. 


INTRODUCTION 

As part of the general program to increase the permissible 
oiDerating temperature of gas-turbine cycles and to improve the life 
of critical parts, further analysis of the liquid-cooled turbine 
blade is being conducted by the NACA to determine the available 
cooling effect at the anticipated hot spot under present conditions 
of design and to evaluate various possible methods of improving 
cooling of the blades. 

Part I of the current series on "Cooling of Gas Turbines" (ref- 
erence l) presents calculations on the application of air cooling 
fins to the turbine rotor with the object of reducing the rim temper- 
ature to peimit better cooling of solid blades; It was found that 
the root of the blade could be cooled about 800® F below the effec- 
tive gas temperature, but relatively little gain was obtained at the 
critical section some distance from the root. In part II (refer- 
ence 2) it was shown that blade life could be extended by rim 
cooling, and that the effective gas temperature could be slightly- 
increased , 


Part III (reference 3) discusses the circulation of liquids 
tlirough hollow passages in the blade and presents a one-dimensional 
analysis of temperature distribution from the blade tip to the rotor 
axis. The effects of variations in blade length, length of cooling 
passage, type of liquid coolant, anci rate of coolant flow were 
studied. The temperatures calculated (one-fifth the gas temperature 
for water cooling) were those appropriate to points near the cooling 
passages. 

The purj)ose of the present detailed study is to investigate the 
temperature distribution encountered in the relatively long, thin 
trailing section between the coolant passage and the trailing edge 
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of uhe liquid-cooled blade. This section, particularly the trailing 
edge, seems most lilcely to be critical. The blade design and 
cooling conditions used are the same as used in reference 3 except 
that only the blade section bounded by the rim, the tip, the wall 
oi* the coolant passage, and the trailing edge is considered. The 
analysis is carried out in several parts, with assumed shapes or 
sections that approximate the trailing portion of a typical blade 
as follows: a one-dimensional analysis of a rectangular section 

through the outer portion away from the root, a similar one- 
dimensional analysis of a wedge-shaped section, and a three- 
dimensional analysis of a right prism for the entire span of the 
blade . 

The following sections of this report present the development 
of the various solutions for temperature distribution in terms of • 
the dimensions, physical properties, boundary temperatures, and 
surface conditions that prevail. In addition, a number of figures 
are presented that illustrate the various effects. of thermal 
conductivity on three-dimensional temperature distribution for 
several values of the effective gas temperature and the combined 
influence of design dimensions, thermal conductivity, and some 
surface conditions on cooling of the trailing edge. 


THEOEETICAL AMLYSIS 
Symbols 

The following symbols, arranged in alphabetical order for the 
convenience of the reader, were used in the calculations: 

a parameter equal to 

b distance from trailing edge to wall of coolant passage, (ft) 

k thermal conductivity of turbine metal, Btu/(hr) (ft) (°F) 

2 distance from blade tip to rim, (ft) 

qj^ heat-transfer coefficient from hot gas to metal, 

Btu/(hr)(sq ft)(°F) 

q^ heat-transfer coefficient from metal to coolant, 

Btu/(hr)(sq ft)(°F) 

t blade thickness, (ft) 
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t mean blade thickness, (ft) 

ti blade thickness at trailing edge, (ft) 

t 2 blade thickness at wall of coolant passage, (ft) 

T temperature of metal, 

Tg effective temperature of hot gas, 

Tj average temperature of liquid coolant, ^ 

Tj, temperature at rim, ^F 

X distance from blade tip to blade element, (ft) 
y distance from trailing edge to blade element, (ft) 
z distance from, median plane of section to blade element, (ft) 

t 

2a angle between sloping sides of wedge 
0 Tg - T, “f 

A linear dimension accented by a prime denotes that the distance 
has been extended by t 2_/2 (for the wedge) or t/2 (for the rectangle). 


Discussion of Simplifying Assumptions 

The turbine wheel with cooling passages is shown, in figure 1 
and a cross section of the blade in figure 2. The approximate blade 
shapes (ri^t prism, wedge, and rectangle) used in this analysis, which 
can be used to approximate any tapered profile, are given in fig- 
ures 3, 4, and 5. 

The equations for the temperature distribution in the blade 
section between the coolant passage and the trailing edge were 
obtained by equating the heat entering and the heat leaving an 
element. In order to simplify the boundary conditions, the cross 
section of the blade under consideration was approximated by geo- 
metric figures having straight sides and areas equal to that of the 
blade, ci’oss section. For the first development, a rectangular shape 
(fig. 3) was used, in which the width was made equal to the distance 
from the trailing edge to the coolant-passage wall and the thickness 
equal to the mean thickness of the blade over this width. For a 
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closei* approximation to the true shape, a wedge shape (fig, 4) was 
used in which the width remained the same as in the rectangle and 
the thickness tapered from the blade thickness at the coolant passage 
to the blade thiclcness at the trailing edge. For both shapes, the 
approximate cooling surface was less than that encountered in tlie 
actual blade because the circular coolant-passage wall was approximated 
by a chord rather than by the true arc. 

The analysis was performed in the following parts: (a) a one- 

dimensional analysis of a rectangular section through the outer 
portion of the blade, away from the blade root; (b) a one -dimensional 
analysis of a wedge-shaped section through the outer portion of the 
blade; and (c) a tliree-dimensional analysis of a right prism for a 
4 -inch blade length. 

The following simplifying assumptions were made for these parts; 

(1) In parts (a) and (b), the heat flow to the rim was negligible 
at a section chosen away from the blade root because the length I 
was large compared to the width b (fig. 5). The validity of this 
assumption is demonstrated by the results obtained in the tliree- 
dimensional analysis. 

(2) The heat gained from the trailing edge can be accounted for 
by assuming the width to be extended by a distance equal to one-half 
the blade thickness at the trailing edge tj^ (fig. 6). The edge 

OQ at temperature T gained some heat. The extended surfaces 00* 
and QQ' were at nearly the same temperature T and no heat entered 
the end 0*Q', therefore these surfaces gained practically the same 
amount of heat as the actual exposed end. In part (c) especially, 
the boundary condition was much simplified. This assumption remained 
valid as long as t/2 was small compared with b. This validity is 
demonstrated for part (a) . In like manner, the blade length was 
extended a distance equal to t/2. The length and the width extended 
by the distance t/2 are denoted by b' and Z'. respectively. 

This assumption was used only in parts (b) and (c ) . 

(3) The temperature of the blade at the rim was constant at rim 
temperature. If this constancy is not maintained, it is shown that 
variations from it have little effect on the results except very 
close to the rim. This assumption was used only in part (c). 

(4) The variation of k was negligible and Tg, Tj, q^, and qg 
are constant over the blade and surface passages. Because these 
blades were assumed to be cooled over the entire length and those of 
reference 3 cooled only to within one-sixteenth inch of the tip, the 
comparable curves differed at the tip. 


6 


MCA EM No. ETBlld 


The analysis of parts (a), (h), and (c) are given in the follow- 
ing sections. 


One -Dimensional Analysis of Eectangular Section through 

Outer Portion of Blade, Away from Blade Boot 

Tlie derivation of the formula for part (a) as shown hy figure 3- 
is obtained from: 

Heat entering element from right end = kt d0/dy 
Heat entering element from sides = 2q,^ dy 6 
Heat leaving element from left end = kt (de/dy + d^o/dy^ dy) 
Therefore, 


kt de/dy dy 0 = kt (d0/dy + d^0/dy^ dy) 


or 


^ _ „2 


dy^ 


= a"-0 


where 


a =M 


2^ 

kt 


The boundary conditions are as follows: 
when y = b, 

kt = (lot (T - Tj) 
If Tq - 0 is substituted for T, 


when y = 0, 


^ (Te - T2 - 0) 


k ^ = q.© 
dy ^ 


A solution of equation (2) is 

0 = A cosh a (y + c) 


( 1 ) 

( 2 ) 

(3) 


(4) 

(5) 

( 6 ) 
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where A and c are integration constants. The quantities e 
and A are evaluated to satisfy the boundary conditions at the 
blade edge and in the cooling passage, respectively. 

From equations (5) and (6), when y = 0, ^ 

icA a sinh a c = qj^ A cosh a € (7) 

therefore from equations (3) and (7) 



when t is not too large. This difference between e and t/2 
varies from 0,4 to 0.5 percent for the values of a considered and 
has an inappreciable effect 0.005 foot or more inside the trailing 
edge . 


When y = b, then from equations (4) and (6) 


Aa sinh a (b+c) = — [Tg-Tj - A cosh a (b + e ) ] 


or 


A = 


T (Te - Tl) 


a sinh a (b + c ) + — cosh a (b + e ) 

k 


(9) 


( 10 ) 


Tlierefore, ^ 


0 = 


9-0 

(Te - T^) cosh a (y + e ) 

a sinh a (b + € ) + — cosh a (b + e ) 

k . 


( 11 ) 


One-Dimensional Analysis of Wedge-Shaped Section 
through Outer Portion of Blade 

It has been demonstrated (reference 4, equation (53)) that a 
solution for the temperature distribution along the y-axis in a 
wedge-shaped segment (fig. 4) may be expressed 

0 = A'Jo(in) + BiHo(in) (12) 

where A' and B are arbitrary constants and Jq and Hq are 
Bessel's functions of two kinds, zero order. 
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The variable |a is defined 


= 40^ 


- , /l - tan a\ 

y h^^- 2Tan"a )_ 


where 


and 




k sin a 


\ _i ^2 - ■tA 

a = tan 1 ; 

The boundary conditions are, when y' = 0 

de 


dy* 


= 0 


when y' = b' 


or when (i = i-ij^ 


and when ij. = iig 


when y' = 0 


when y ' = b ' 


de 


k ~ = To (Te - Tz - 0) 


d.\x 


= 0 


= (Te-T, -0) 


hX = 2C \ 


Jti (l - tan a) 
y 2 tan a 




d6 

Equation (13) can be differentiated and put in the form — , 

de _ ui de 
dn " 2C^ dy’ 


Therefore when 


then 


de 

dy’ 


= 0 


By virtue of the properties of Bessel’s functions, 

dJo (y)/dy = “Ji (y) 


(13) 

(14) 

(15) 

(16) 

(17) 

(18) 

(19) 

giving 

(20) 


( 21 ) 
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and 


therefore 


and when 


then 


dn 


dH^ ( 7 )/dy = -% (y) 

= -A'iJi (i^I) + B Hi (in) 


1^ = 0 

d|i 


(22) 


(23) 


® - Hi (ini) 

In order that eq.uation (17) Diay he satisfied 
2 

k ^ [-A'i Ji(in2) + BHi(in2)] = aoC(Te-Tz) - A'Jo(in2) - Bi Ho(in2) ] 

(24) 

In standard form eq.uations (23) and (24) can he written, respectively 
[iJi (ini)] A' - [Hi (ini)] ® = ° 






B = Te - 


The values of A' ajid B can easily he found and suhstituted into 
equation (12), giving 6 in the form 

qoM-2 

(Te-T^) [Hi(im)Jo(in) + iJi(im)iHo(in)] ' 
e = — __ (25) 


[iJl(ini)Hi(in2)] - [Hi(ini)iJi(in2)3 + X 

SC'^k 


where X is 


[Hi(lni)Jo(in2) iJl(it^l)iHQ(in2)J 


Three-Dimensional Analysis of Right Prism for 
4 -Inch Blade Length 


The derivation of the formula for part (c) as shoi-m in fig- 
ure 5 is obtained from ^ 
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Heat entering element from top = -k djmz 


St 


Heat entering element from right end = -k dxdz 

St 

Heat entering element from front = -k dxdy j- 


dT 

55^ 


Heat leaving element fi-om hottom = -k dydz 1 


St . S^T 3 . 


_ / ^ip 

Jleat leaving element from left end = -k dxdz | ^ + g— g- dy 

/St S^T \ 

Heat leaving element from hack = -k dxdy + -— 5 - dzj 

\oz 32^ j 


When the heat entering is equated to the heat leaving 


St St St /^St S2t '' 

-k dydz ^ -k dxdz -k dxdy = -k dydz i ^ ^ dx 

t> ^ OX j 


, , , /St S^t , \ , /ST s<^T , \ 

-k dxdz i^gj: + dyj -k dxdy ^ azj (26) 


/St S^t 


or 


0 = ^ ^ ^ 

Sx2 Sy2 Sz2 


(27) 


The boundary conditions can he more simply handled if the gas 
temperature is taken as a reference temperature rather than the 
metal temperature; that is, if the suhstitution is made 


a = Tg - T 


(28) 


or 


T = Tg - a (29) 

Equation (27) then becomes 

S2e S2e S2e 

0 = ^ + — -p + — J (oO) 

Sx'"^ Sy'“^ Sz^ 

The origin of the coordinates chosen is jjidicated in figure 5. 
The axis z = 0 is located in the median plane of the right prism; 
from considerations of sj^imaetry there can he no heat flow across the 
median plane. 
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The boiindary conditions that must he satisfied are 


when x' = 0 


when y' = 0 


when z = 0 




(31) 


(32) 



(33) 


when x' = 


when y* = b* 


when z 


t 

2 


0 = Tp - Tv 


^ = % (Tg - Tj - 9) 


k^0 

5T = -'ll ® 


(34) 


(35) 


(36) 


A convenient form for the solution of equation (30) is 

CX CD 
1 1 

X (Ajjjn cos %x' cosh t^^any' cos PmZ +Bnm cosh EuinX' cos Any' cos Pmz) 


(57) 

where A, M, N, E, K, and P are arbitrary constants. The rela- 
tion between the constants used as coefficients for the variables is 

^'Wi^ = Pm^ + (38) 

Pmn^ = Pm^ + ^n^ (39) 


where n = 12, 3 4, ..., and m = 1, 2, 
tion oi equation (o7), the boundary condit 
tions (31) to (33) are satisfied. 


3, 4, .... By the selec- 
ions expressed in equa- 
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and Bmn are derived. These 


1 

sinh ^^h’ + q,o cosh 


t 

2 1 1 

^ ~ h{b' , . Bin p t ^ t cosh E^jn 1' 

1 + — ?rr " r -; -*^111 ^ Bin ^ 

1 + ^ 

2P — 

Values of n froia 1 to 14 vrere inserted and the temperatures cal- 
culated from equation (37), 

In appendix B it is shoTO that sufficiently accurate results 
(within 1°) can he obtained by using only the first value of m; 
that is, m = 1 and dropping the terms involving cos P 2 z, 
cos Pj z, etc. 


In appendix A, values of A ^n 
values are 

. V 

(_l)n-l am ^ (Tg _ 

^ - «(2n-l) _ t . 0-0 

Pm g- sin 2Pm ^ 

2 "° — 


sin X.nt‘ 4 (Tg - Tr) 


m 


APPLICATION OF ANALYSIS 

The following assimiptions were made in applying the previous 
results to specific numerical calculations: 

1. The gas flow was 55 pounds per second. The heat-transfer 
coefficient corresponding to this gas flow qp was found to be 
222 Btu per hour per square foot per °F. 

2. The liquid coolant, water, had an average temperature of 
200^ F and a flow rate of 6.42 pounds per second for the entire rotor. 
The heat-transfer coefficient corresponding to this coolant flow q^ 

was calculated (reference 5) as 2370 Btu per hour per square foot per °F. 

3. Average themial conductivities k of 15, 40, 60, 80, 100, 

120, and 210 Btu per hour per square foot per °F were used. 

4. The turbine had 55 blades. 

5. The blade was liquid cooled over its entire length of 4 inches, 
by two cooling passages 0.25 inch in diameter. The blade chord was 
1.188 inches long. 
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6. Some numerical values assumed were 
b = 0.05C0 foot 

b' = 0.0550 foot (rectangular section) 

b' = 0.0515 foot (wedge-shaped section) 

I = 0.3333 foot . 

2* = 0.3383 foot 

t = 0.010 foot (rectangular section) 

t]_ = 0.003 foot (wedge-shaped section) 

t 2 = 0.021 foot (wedge-shaped section) 

Tj. = 284° F (all values of thermal conductivity) when 
Tg = 2000° F 

= 331° F (all values of thermal conductivity) when 
Te .= 3000° F 

Tp = 426° F (all values of thermal conductivity) when 
Te = 5000° F 

f 

The theoretical temperature distribution in the trailing segment 
of the turbine blade was determined by a three-dimensional analysis 
using equation (62) in which the blade cross section was approximated 
by a rectangle. The temperature distribution was found in two planes 
representing the maximum and the minimum temperatures for the z 
axis; the first plane was located at the side of the rectangle and 
the second was on the median plane througli the section. The tempera- 
ture was plotted eigainst radial distance from the wheel axis for 
intervals from the trailing edge to the wall of the coolant passage 
(fig. 7). At the wall, the series converges too slowly to be useful, 
therefore the curve was based on a one -dimensional calcu3.ation (equa- 
tion (ll)) and the behavior of the other curves of the family. Sim- 
ilar distributions were calculated and plotted using theimal conduc- 
tivities 120 £md 210 Btu per hour per foot per °F (fig. 8). 

The next step in the analysis was to determine the effect of using 
a rectangular cross section as an approximation of the true blade 
shape. The complexity of a solution using the true blade cross 
section eliminates an exact comparison. A wedge-shaped cross section 
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represents the closest approximation with planes as boundaries and 
this shape was used to indicate the trend. A three-dimensional 
analysis of a wedge-shaped cross section was beyond the scope of 
this report and quite unnecessary because on the outer three-quarters 
of the blade, the metal temperature does not change with the radius 
and the comparison was made on a one-dimensional basis for tempera- 
tures along the center line of a section through the outer portion 
of the blade. The comparison is presented as a plot of metal tem- 
perature against distance froia the wall of the coolant passage for 
an equivalent rectangular and wedge-shaped section (fig. 9). 

The thermal conductivity of the metals used in blade construc- 
tion has considerable influence on the temperature distribution; the 
high-tensile-strength alloys currently used have a relatively low 
thermal conductivity. The effect of thermal conductivity was 
investigated using a one -dimensional analysis of a section through 
the outer portion of the blade for average thermal conductivities 
ranging from 15 (steel) to 210 Btu per hour per foot per °F (copper). 
Although the upper range of conductivities is impracticable for 
actual use in blade construction, it does indicate the trend and the 
magnitude of temperature changes with changes in thermal conductivity. 
The comparison is presented as plots of metal temperature against 
distance from the wall of the coolant passage for. the rectangular 
siiape (fig. 10(a)) and the wedge shape (fig. 10(b)). The effect of 
a change in thermal conductivity on trailing-edge temperature was 
determined comparatively using the same values of conductivity for 
equivalent rectangular and wedge shapes. One station on the section 
center line at the trailing edge was used for each thermal conduc- 
tivity (fig. 11). 


RESULTS AND DISCUSSION 

The theoretical three-dimensional temperature distribution for 
a turbine blade of rectangular cross section is presented for effec- 
tive gas temperatures of 2000^"’, 3000°, and 5000° F in figure 7. The 
temperature falls off in increasing increments as the distance to 
the cooling liquid decreases. At a constant distance from the wall 
of the coolant passage, the temperature remains substantially constant 
for approximately tliree-quarters of the blade length. The tempera- 
ture beyond this point falls off rapidly and reaches the rim temper- 
ature at the base. The average temperature at the trailing edge of 
the blade was reduced 7.14 percent below the gas temperature at 
2000° F and 7.6 percent below the gas temperature at 5000°. The 
temperature differential between the side and the center of the 
blade (z = t/2) was 3.6 percent of 6, negligible at the trailing 
edge, and at the maximum point did not exceed 5 percent of the metal 
temperature . 
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The temperature distrj.hution at a gas temperature of 2000° F 
for rectangular blades having thermal conductivities of 120 and 
210 Btu per hour per foot per °F is presented in figure 8. With 
increasing thermal conductivity of the blade metal^ the temperature 
tends to decrease at the trailing edge and increase near the coolant 
passage. The temperature at the trailing edge is approximately 
32 percent below the gas temperature at a thermal conductivity of 
120 Btu per hour per foot per °F and 36 percent below the gas tem- 
perature at a thenaal conductivity of 210 Btu per hour per foot 
per Op. With increased thermal conductivity, the influence of rim 
cooling is more effective and- the radial distribution lines fall 
off more rapidly. Tlie temperature differential between the sides and 
the center of the blade decreases with increased thermal conductivity. 
For the higher conductivities, this differential is too small to be 
sho\m in the figures, inasmuch as it ranges from 2° to 7° F. 

An investigation of a rectangular cross-section blade indicates 
the temperature distribution in a blade having the cross section 
shown in figure 2 but errors in the temperature distribution 

®^ist. The trend will be indicated, however, if the rectan- 
gular cross section is compared with a closer approximation to 
the true shape; namely, a wedge-shaped section. Figure 9 presents 
plots of the one-dimensional temperature distribution along the 
center line at the blade tip for a rectangular and a wedge section 
at a gas temperature of 2000° F and a thermal conductivity of 15 Btu 
per hour per foot per °F. The temperature distribution for the 
wedge shape has a slightly steeper slope than tliat for the rectan- 
gular section and the temperatures at the trailing edge and the wall 
of the coolant passage are lower. (Part of this lowering is due to 
the additional thickness of the cooling surface.) 

The effect of the thermal conductivity of the blade metal on 
one-dimensional temperature distribution is shown in figures 10 
and 11. Figure 10(a) presents the temperature distribution along 
the center line of a rectangular section through the outer portion 
of the blade for an effective gas temperature of 2000° F and various 
thermal conductivities. At low thermal conductivities, the temper- 
atui-e falls off sharply from the trailing edge to the coolant 
passage. With increasing thermal conductivity, the temperature 
decreases at the trailing edge and increases at the wall of the 
coolant passage, which in general tends to equalize the temperatures 
across the width of the blade. The reduction of the temperature 
at the trailing edge increases with an increase in the thermal 
conductivity. Figure 10(b) shows the temperature distribution 
j. or a wedge section for the same conditions imposed on the rectan- 
gular section. The general form of the curve is similar to that 
for the rectangular section; hovrever, the temperatures at the 
trailing edge and at the wall of the coolant passage are lower and 
the drop in temperature for a given increase in conductivity is greater. 
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The variation in temperature with an increase in thermal con- 
ductivity for rectangular- and wedge-shaped sections is shown in 
figure 11. The point selected for comparison was located on the 
center line of the section at the trailing edge. The curves show 
that the blade temperature at this point may be considerably lowered 
by using metals having a lii^ thermal conductivity and that for the 
same conductivity, the cooling is substantially greater for the 
more representative wedge-shaped section. This effect is partly due 
to the loss of effective area at the cooling passage where the thick- 
ness of the equivalent rectangular section is considerably less than 
that of the wedge section. At high values of thermal conductivity 
a representative teraperature distribution can be obtained from the 
rectangular section by assuming the mean thickness t equal to the 
value of t 2 for the comparative wedge. The total cross-sectional 
area then is no longer equivalent. The same care in selecting 
appropriate dimensions extends the useful range of the nondimens ional 
chart (fig. 12). . 

The nondimens ional cooling ^ at the trailing edge for 

l0 - 

the rectangle is shovn as a function of ah’ and in fig- 

ure 12. Equation (ll) was used with cosh ac ^1* It is seen from 
figure 12 that for ah* ^ 3, the cooling is less than 10 percent 
for the entire range of qo/ak. In order to obtain 50 percent 
cooling even for quite large values of qo/ak, ab ' must be less 
than 1. Points for the extreme values of k, 15 and 210 Btu per 
hour per foot per °F, are shown on the figure. Theoretically ab' 
could be reduced by decreasing q^^ and b' or by increasing k 
and t. Under current conditions, not much can be done about qj_ 
and k. The curves then indicate short thick trailing edges instead 
of long thin ones for good cooling. From the curves of figure 12 
it can be deduced that with materials of low thermal conductivity 
an infinite increase in the coolant-passage heat-transfer coeffi- 
cient relative to that obtained with water will have small effect 
on the trailing-edge temperature. 


CONCLUSIOIK 

The follovring general conclusions can be Arami from the analysis 
of the cooling characteristics of the trailing part of a liquid- 
cooled turbine blade: 

1. For blade materials having low thermal conductivity, the 
computed three-dimensional temperature distribution of the equiva- 
lent rectangular trailing section indicates that the rim-cooling 
effect is significant for only a short distance from the blade root. 
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the trailing-edge temperature closely approaches the effective gas 
temperature^ and the temperature gradient normal to the surface 
exposed to the combustion gas can generally be neglected. 

2. The three-dimensional and one-dimensional solutions of the 
equivalent rectangular trailing portion give almost identical tem- 
perature distribution in the blade section away from 'the influence 
of rim'coolingj the one-dimensional trealanent of the wedge or 
i'ectangular section is therefore sufficient for most investigations. 

3. For blade materials of very low thermal conductivity, the 
one-dimensional temperature distribution computed with the equiva- 
lent rectangular section is only sliglitly different from that of 
the more representative wedge section. 

4. The equations for the' one-dimensional rectangular solution 
permit construction of a simple chart, which gives the value of the 
trailing-edge nondimensional cooling parameter with any possible 
combination of the design and operating variables considered in the 
analysis, 

5. For the blade shapes considered in this analysis, and with 
materials of low thermal conductivity, an infinite increase in the 
coolant-passage heat-transfer coefficient relative to that obtained 
with water will have small effect on the trailing-edge temperature. 

6. For the blade shapes considered in this analysis, trailing- 
edge temperature may be considerably reduced with a material of high 
thermal conductivity and to a lesser extent by locating the cooling 
passage closer to the trailing edge. 

7 Large reduction in trailing-edge temperature results when 
the blade shape is altered to provide a short thick trailing section. 


Aircraft Engine Eesearch Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio. 
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For thes9 yalues of N equation (40) assumes the form 
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sin Afjb* 
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f eln2Anb' \ 

V 2\ih' y 


( 61 ) 
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APPENDIX B 


EVALUATION OF THE TEEMS v/BERE m > 1 

Hie temperature distribution in the turbine blade under consid- 
eration is determined neglecting the smaller terms when m > 1 and 
using equation (37) in the form 

PC; 

9 = cos P]^2i^n (A^jj cos N^x* cosh +3 -|-n cosh E]j^x’ cos Any') 

^ (62) 

This equation will be valid for all practical purposes when the 
thickness of the blade is small compared with the width and the 
length. It is necessary^ in deaionstrating this statement, to expand 
the double series given by equation (37) and inspect the terms con- 
taining the function cos PmZ. The expansion of equation (57) with 
the subscripts in the order mn assumes the form 


0 = (All cosh Mp;j^y' cos PpZ+A2]_ cosh M2]_y' cos P2Z + ..,) cos NpX* 

+(A]_2 cosh Mp2y' cos PpZ + A22 cosh M22y' cos P2Z + •••) cos N2X' 

+. . . + (A]_j^ cosh M-|riy ' cos P2_z +A2n cosh M^ ^ y ' cos P2Z+ ..,) cos Nj^x‘ 
+(Bii cosh Ej^px’ cos Ppz +B2X cosh E21X’ cos P2Z + .,.) cos Apy’ 

+(Bp2 cosh Ep2x' cos Ppz+B22 cosh E22X' cos P2Z+ ♦..) cos \ 2 J' 

+ ... + (Bpn cosh EpnX' cos Ppz +B2n cosh E2nX' cos Pn^ + .*.) cos XnX' 

( 63 ) 

If all the other terms in the series are arranged so they become 
coefficients of the function cos P^z, 
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/^ A A21 cosh M217' \ 

e = All coBh Miiy' UoB PiA . All cosh Miiy' ^2' * • ' •] 


/ A' 

+ Ai 2 cosh Mi2y' (cos P^z + j 


A22 cosh M22y' 


\ 


12 cosh Mi2y’ 


7 cos P2Z + ,.,J cos W2X* + . . . 


A -u X. , -r. A2n cosh M2ny' ' 

+ Ain Mljiy' Pl^ + Am cosh Mmy' ^2" * • • •} 

■r, X -rx . ' ^ ®21 cosh E 2 IX’ \ 

+ Bii cosh Riix« ;^C0S PiZ + g- ■ cos P2Z + . . .j cos \iy' 


+ Bj^ 2 cosh 


•n t / TA ®22 cosh R22X* 'i 

E12X ^cos Piz + B12 ' cosh Ri2X» P22 + • • •] cos ^gy' 

T> A. TA . AA ®2n cosh P2nX' \ 

A Bm OOBh EmX' I COB PiB . g --- COB PjZ . . . . j OOB V' 

( 64 ) 


From equation ( 64 ) it can “bo seen that if the coefficient of the 
function cos P2Z in any series is negligible compared vith unity^ 
the terms subsequent to cos Ppz may be disregarded for practical 
application of the equation. Because the series is convergent for 
ascending values of n^ it is only necessary to inspect the coef- 
ficient of COB Pgz for a value of n = 1. 


V/hen the coefficient of cos P2Z in the Ajjjq series is equated 

to 0 and the values of Ajj^ determined by equation (59) are substi- 
tuted^ 


_ t‘ / t\ 

sin P2 2 / 2 Pp 

- .1 + ^ 


0= 


1*2 5- 


\ 


2Pi - 


„ t / 
sin Pi g- / sin 2P2 g 
il + 


2 / (lsi 4 ii sinh Mub’ + qo cosh Muh’) 

t \ (kJ'121 sinh M2ih’ +qQ cosh M2ih') 


j. ' J 

V 


®2 I J 


cosh M2 i 7* 
cosh Mppy' 


( 65 ) 



26 


MCA EM No. ETBlld 


or 




sin 'So% f sin 2Pi 

1 ( 1+ £■ ! 

V 1 \ 9V 1 , 

■*^2 2 \ ^•*^1 2 / (isi'ill taiih Mnii' +q.o) cosh Muh’ cosh M 2 U‘ 


sin 1*1 2 / 


(1cM 21 tanh M 2 ih' +<lo) cosh M 2 ]h' cosh Mi i.y' 


■1 2 


X66) 


For the "blade under consideration. 


h‘ = 0,055 ft 
k = 15 Btu/(hr)(ft)(°F) 

1* = 0.3383 ft 

Qj^ = 222 Btu/(hr)(oq ft)('^F) 

= 2370 Btu/(hr)(sq ft)(°F) 
t = 0.01 ft 
From equation (50) 


= 53.74 

P 2 = 632.99 
From equation (43) 

Ni = 4.64 ft-1 
From equation (38) 

Mil = 53.94 ff^ 

M 21 = 652.99 ffl 

The following circular and hyperbolic functions apply for these 
values : 
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Sin Pq 1 = 0.2654 

cosh Mj^j^b 

sin Po i = .0234 

2 2 

cosh Mgqb 

sin 2P^^ 1 = 0.5119 

cosh Mqqy 

sin 2Pg 1 = .0467 

cosh Mgqy 

tahh Mqqb' = 0.995 

cosh Mqj^y 

tanh M 23 _b* = 1 

cosh M 2 qy 


= 9.7425 
= 6.5907 X 10^^ 

= 1 when y' = 0 
= 1 when y’ = 0 
= 7.4513 when y' = .05 
= 2.7822 X 10^^ when y' = ,05 


If the given values in equation (66) are substituted, the following 
coefficients are obtained: 


y' = 0 0 = -5.738 X lO-l”^ 

y' = 0.05 0 = -2.142 X 10’^ 

y‘ = b' 0 = -3.882 X 10"^ 

The values of 0 obtained for the blade under consideration show 
that the coefficient of the cos Pgz term has a maximum value at 
y’ = b’ and is then equal to 0.4 percent of the first term and as 
such may be disregarded for practical applications. 


The effect of the second term of the series can be inves- 
tigated by equating the coefficient of the cos P2Z term to x and 
the previous process repeated. 


E 


t / t\ 

Bin P2 o- / sin 2P-i 5-'. 

r^l 

^2 2 \ ^^1 2 / cosh E]^pZ* cosh ^21^* 


^1 2 /. 

^l| 


1 + 


t\ cosh Bp-iZ' cosh BtiX' 
sin 2?2 p \ 


(67) 


2P2 2 / ' 


The application of equation (67) to the blade under consideration 
requires the following values: 


From equation (48) 


\l = 25.64 ft-1 
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From equation (39) 


Ell = 59.54 ff ^ 

E 21 = 633.51 ft~ ^ 


cosh EjjjjjX' 
■becomes 


With the values of S, 

may be approximated by 


‘mn used, the cosh U nm l' and 

e^mni ’ 


and equation (67) 


t / 

sin Pg 2 i 2 Pp 2 

■ (1 + ^ j 


r = 


p - 

■2 2 


2P ~ 

e(^ ' )(%1-I^2l) 


sin Pi - / sin 2P2 

— — h + 2 i 

1 \ 2Po ^ / 

2 \ ^^2 2 y 


( 68 ) 


Pi 


When the given values in equation (68) are substituted for a 
range of values of the following values of V are found: 


X' 

(ft) . 

Distance from rim | f 

(«) 

0.005 

0.3335 j 

.3283 

.01 -0.4643 X lO'’; 

.3333 

.005 ! - .8227 X 10"2 

.3383 

0 1 - .1451 X 10-P 


The coefficient of the cos term increases as the rim is 

approached and reaches a value of 1.5 percent of the first term at 
this point. For purposes of this report^ the temperature will not 
he calculated withn.n 0.005 foot of the rim and the cooling liquid. 
With this limitation^ the values of the second and subsequent tenas 
are negligible . 
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Fig. 7a 



Figure 7. - Three-dimensional temperature distribution in rear part of turbine blade 
at thermal conductivity k of 15 Btu/( hr) ( ft) (°F) . (For facility in reading, fig. 7 
(a) is plotted on a different ordinate scale than figs. 7(b) and (c).) 
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0 I 2 3 4 

Blade length, in. 


. (b) Effective gas temperature, Tg, 3000° F. 


Figure 7, - Continued, Three-dimensional temperature distribution in rear part of tur 
bine blade at thermal conductivity k ot 15 Btu/( hr) ( ft)(°F). 
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Figure 7. - Concluded. Three-dimensional temperature distribution in rear part of tur 
bine blade at thermal conductivity k of 15 8tu/{ hr) ( ft) (°F) . 
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(a) Thermal conductivity, k, 120 Btu/( hr)( ft) (°F) . 


Figure 8. - Three-dimensional temperature- distr i but ion in rear part of turbine blade at ef- 
fective gas temperature Tg of 2000^ F. 
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Fig. 8 b 
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(b) Thermal conductivity, k, 210 Btu/(hr )(ft) (®F) • 


Figure 8. - Concluded. Three-dimensional temperature distribution in rear part of turbine 
blade at effective gas temperature of 2000® F. 
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Figure 9. - One-dimensional analysis of the effect of shape on prevelanet blade tem- 
perature. Effective gas temperature, Tg, 2000° F; thermal conductivity k, 15 Btu/ 
lhrHft)(°F). , . , 
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Fig. I Oa 
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(a) Rectangular section. 


Figure 10. - One-dimensional solution for pervalent blade temperature at effective 

gas temperature of 2000® F. 
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(b) Wedge-shaped section. 

Figure 10. - Concluded. One-dimensional solution for prevalent blade temperature at 

effective gas temperature Tg of 2000° F. 
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Figure II. - Temperatures of trailing edge obtained from the one-dimensional solutions 

Effective gas tem^rature, 2000P F. 









